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I. INTRODUCTION 


In this paper, following a suggestion of Zagier Q, we generalize the “number theoretic” partition function, whose 
statistical mechanical properties have been studied by Knauf iiBi, by introducing a new parameter x £ Rj. 
Both the “canonical” and “grand canonical” partition functions of Knauf arise, for certain values of x. More generally, 
as explained in section im these models may be regarded as one-dimensional spin chains of length k , and the new 
parameter x allows us to derive recursion relations on the length of the spin chain. These relations are simple 
generalizations of the Lewis three-term equation that has been extensively studied in number theory 00 Next, 
section ITTT1 explores some simple consequences of spin orientation (“spin flip”) transformations for partition functions 
and expectation values. These results, along with the recursion relations, are our main tools. They are used, in 
various ways, in the succeeding sections. In section ITvl the recursion relations are shown to imply a very simple 
and direct connection between the transfer operator studied by Prellberg ns ns in a model of dynamical systems 
with intermittency and the generalized partition function. This is one of our main results. We examine some of its 
consequences. In particular, we prove that all models have the same free energy (and hence the same thermodynamics), 
independent of the value of x, and show that the known spectrum of the transfer operator implies that the correlation 
length satisfies the prediction of scaling theory. In sections El and eh we use the recursion relations and spin-flip 
behavior to calculate certain spin expectation values for both finite and infinite spin chains. Some new and interesting 
features arise. In particular, the expectation value of certain spin clusters and an independent spin at arbitrary 
distance is shown to be independent of the direction of the spin, at all temperatures above the transition. Thus, in 
this sense, the spin cluster removes the spin asymmetry of the system. This behavior appears to be related to the 
mechanism of the phase transition. 

All of our results are rigorous. Since they concern certain weighted averages over Farey fractions and relate to 
the Lewis equation, they may be of interest to mathematicians. Therefore we have included a few explanations and 
definitions in an attempt to make the paper more accessible to those unfamiliar with statistical mechanics. 


II. DEFINITION OF THE PARTITION FUNCTION 


In this section we define the generalized partition function, and show that it satisfies a recursion relation. This 
relation is one of our main tools for proving new results. 

Let the matrix M k be any product of k matrices Aq := (} j) and A\ := (* *), 


M k = 



We can regard any such product as a one-dimensional chain of length k. If the itli matrix (1 < i < k) is A 0 we identify 
it as a spin pointing up at the ith site in the chain, and likewise as a spin pointing down if it is A±. Therefore each 
matrix M k corresponds to a definite configuration of the k spins. 

Next we extend the Knauf model Q (see also EH and 0 ) by introducing a family of partition functions parametrized 
by the variable x > 0 [H 


Z k (x,P) := ^{cx + d) 2/3 , 


( 1 ) 


where the sum runs over all 2 k permutations of the product of the k matrices Ao, A\. 

Next, setting M k+ 1 = M k Ao or M k+ 1 = M k Ai gives rise to a chain of length k + 1. Now since M k A^ 

and M k Ai = ( “ ), we find from Q the recursion relation 



Zk+i(x,0) = (1 + x)~^Z k ^JL-^ + z k {x + 1,/3) (2) 

with the initial condition Zq{x, (3) = 1 (i.e. Mq = (JJ)). The variable x £ Rj is a parameter which changes the energy 
(E k = 21n(c:r + d )) of each spin configuration (i.e. each matrix product M k ) and /3 £ R(| is the inverse temperature. 
However, the thermodynamics is independent of x , as we will see. 

It is convenient to define, as in number theory E3 , the action of the matrix M = on any function f(x) 

f{x)\M := (cx + d)~ 2l3 f ■ (3) 
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For example, consider the action of the matrix Aq on the constant function 


l{x)\A 0 = (1 +x) 2/3 , 


( 4 ) 


where l(ir) = 1. 

It is easy to check that our partition function Z k (x,(3) can be written as 

2 k 

Z k (x,P) = ^2l(x)\Mi, ( 5 ) 

2=1 

where M,; = TLj =1 A Tj ( i ' ) with Tj(i) € {0,1}. Note that each defines fractions f and ^ at level k of the Stern-Brocot 
tree 0 ; thus the level corresponds to the length of the spin chain. The subset of these fractions between zero and 
one are called Farey fractions. They are generated by the products which start with Aq 0 - 
In the following, we make extensive use of an abbreviated form of © 

Zk{x,P) = 1 (x)\(Aq + Ai) k = 1 {x)\Aq(AoA i) k 1 + l(a;)|j4i(j4o + A\) k 1 , (6) 

where the addition must be applied after the multiplication of the matrices! 

The Knauf “canonical” partition function Zj?(s) (see jlj for the definition and note that s = 2/3) is equal to 

Zff = 1 ( x )\Ao(Aq + Ai) k \ x =o. ( 7 ) 

Similarly, the “grand canonical” partition function of Knauf ;2] corresponds to 0 with x = 1 on the right hand side. 
Let 


Zk(x,/3) l(x)\A 0 (A 0 + Ai) k . (8) 

0 is then a direct generalization of Z^. It is easy to verify that 0 satisfies the recursion relation ©. Using 
\{x)\A\ = l(:r) and 0 we get 

k -1 

Z k (x,P) = Z k _ 1 (x,P) + l(x)\A 1 (A 0 + A 1 ) k - 1 =l + J2 z i(x,P). (9) 

i=0 

Thus 0 relates two partition functions satisfying the recursion formula 0 with initial conditions Zq(x, (3) = ( l+x )~ 2/3 
and Zq{x,[ 3) = l(x). 

In the following, we use Z k {x,P) exclusively, because of its direct relation to previously studied spin chains. It is 
possible to obtain similar results for Z k (x,/3) as well. We show in section llVl that all of the partition functions 0 
have the same free energy, i.e. their thermodynamics is independent of x (for the definition of free energy, cf. 13611 1. 
For use below, we note that it is straightforward to express Z k {x,f3) similarly to Z k {x,/3) in 0 by 

Z k (x,P) =^2((a + c)x+(b + d))~ 20 , (10) 

where a,b,c,d now correspond to the matrix elements of M k -i- Note, however, that (see 0 ) M k is now always of 
the form AoM k -i■ Thus, setting x = 0 (the “canonical” case), one sums over all Farey denominators, given here by 
b + d, at level k. Letting x = 1 to obtain the “grand canonical” partition function thus corresponds to summing only 
over the “new” denominators (a + c + b + d) at the next level. Now the Farey fractions at each level k are composed of 
“old” fractions that arose at lower levels and “new” ones from level k. Thus the “grand canonical” partition function 
at level k can be written as a sum over “canonical” partition functions at lower levels, i. e. over all “canonical” chains 
of shorter length. This is the opposite of the usual situation in statistical mechanics, and the reason why we put the 
names in quotes. 


III. SPIN ORIENTATION TRANSFORMATIONS 

In this section we consider the consequences of the spin-flip transformation generated by the matrix P = (^ J). 
Specifically, we investigate the effects of P on the partition function Z k (x,/3) and some related functions useful in 
calculation expectation values. 
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The action of P on a function f(x) is 


f(x)\P = x 2,3 f{l/x). 


( 11 ) 


Note that the matrix P simply exchanges the spin orientation, e.g. the matrix Aq and the matrix A\ satisy 


A 1 =PA 0 P. ( 12 ) 

Since P 2 = 1, iTHU) in fact implies that A 0 and A 1 are conjugate. Note that a function f(x) satisfying ITT1) (i.e. f(x) = 
x 2l3 f{ 1/a;)) can be called even, since, using the substitution e y = x (recall that x > 0 herein) to define g(y) = e /3 /(e y ), 
<tTTI becomes g(y) = g{—y)- 

Now our initial condition Zq(x,P) = (1 + x)~ 2 ^ is easily seen to be even. Consequently, for all k > 1, x £ R + and 
/3 £ Rq" the partition function Z k (x) is even 

Z k (x)\P = (1 + x)~ 20 \{A Q + A{) k P = (1 + x)- 2/3 \P 2 (Ao + Ai) k P = Z k (x). ( 13 ) 

In the last equality we used the evenness of our initial condition and the fact that the set of all terms in (Aq + Ai) k 
is the same as the set P{A 0 + Ai) k P. Note also that 10 implies that the partition function Z k (x,p) is “almost” even 
for [3 < P c and k —> oo, since both Z k (x,P) — 1 and Z k (x,P) diverge in this limit, and the former is even. Finally, 
(fTT]) and m show that Z k { 1), the “grand canonical” partition function, is actually invariant under the spin-flip 
transformation. This corresponds to the absence of odd-spin interactions in this model, as will be discussed below. 

Now consider the terms in 0 . Using the evenness of our partition function we can write Z k -i(x) |^4o = Z k -i(x)\P Aq 
and Z k _i(x)\Ai = Z k _i(x)\P A\. Thus 

(i + x)~ 2l3 z k _ 1 = x ~ 20 z k _i [y-y-iP'j ( 14 ) 

and 

Z k -i(x + l,/3) = (1 + x) 2 ^Z k -i ^-j——j—^, P^j (15) 


for all k > 1, x £ R + and P £ Rg . Combining (thill . (fTTHl and 0 (which, as mentioned, also holds for Z k (x,P)) gives 
us four different possible recursion formulas. For instance 


Z k ( x) 


C x + 1)~ 20 


Z k -i 



+ -Zfc-i 


x + 1 


(16) 


which we will use in section UYI 

In addition we can see that the matrix P can be put in front of any matrix Aq or Ai in the expression (l+a:)^ 2 ' 3 1 (^4o + 
Ai) k without changing the partition function Z k (x) (for example (1 + x)~ 2l3 \{A Q + Ai) l (P A 0 + Ai){A 0 + Ai) r = 
(1 + x)~ 2 P\(Aq + A\) k for any k,l,r > 0 such that l + r + 1 = k). On the other hand if we put the matrix P after any 
matrix Aq or Aj we get a new function. Let 


Z k r {x) — ^(1 + x ) 2/3 |(Alo + Ai) l {A 0 + A 1 P)(A 0 + Ai) r 

and 

Z' k T {x) = 2 ^ x ) 2/3 |( j4 o + Ai) 1 (AqP + Ai)(Aq + Ai ) r , 
with l + r + 1 = k. Using d we then have 

Z^(x) = Zi(x)\A 0 (A 0 + A 1 y 

and 


Z l k i ^x) = Z l {x)\A 1 {A 0 + A 1 ) r . 


(17) 


(18) 


(19) 

( 20 ) 


The arrows | and | refer to the interpretation of Aq and A±, as up and down spins, respectively. Thus m and 
motivate the notation in and m- In addition note that 


Z^[x) + Z^{x) = Z k [x). 


( 21 ) 
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We will make use of these functions to calculate expectation values in section El 

Note that, as is often done in statistical mechanics, if we fix the spin at one position (or spins at several positions) 
and sum over the rest, as in m or PHi . and then divide by the partition function, the result is an expectation value, 
since the ratio is the sum of the probabilities of all configurations with this spin (or these spins) fixed in the way 
specified. This follows because each term in the partition function is the unnormalized probability of the corresponding 
spin configuration. In statistical mechanics, expectation values involving more than one spin are sometimes referred 
to as “correlations” or “correlation functions”, especially when one focuses on their dependence on the distance(s) 
between the spins. 

We conclude with an observation which follows immediately from El and El. The probability of a spin up at 
position l + 1 from the left is equal to the probability of a spin down at the same position for a model with different 
x. It is easy to see that 

Z‘ h U (x) = Z^(x)\P = x~^Z^( l/x) = Z\} r (l/x) 

Z k (x) Z k (x) Z k (x)\P Z k ( l/x) 

for all k > 1, x £ R + and f3 £ Rj. Note that for x = 1 these probabilities are equal. For other values of x , since the 
magnetization (which is essentially the probability of spin up minus the probability of spin down) is zero, the up and 
down spins probabilities become equal when l and r are sent to infinity (see section IVll. 


IV. CONNECTION TO THE TRANSFER OPERATOR 


In this section we demonstrate a direct and simple connection between the partition function Z k (x, (3) and a transfer 
operator for a model of intermittency in dynamical systems associated with the Farey fractions. This connection had 
already been noticed in ijj, but in a less direct setting. Our new result allows us to prove that the free energy (cf. 1331) 1. 
which is given by the largest eigenvalue of the operator, is independent of x, and draw other conclusions as well. In 
particular, the spectrum of this operator has been determined by Prellberg Q, and it follows from his results that 
there is a second-order phase transition for all x. 

To begin, consider the Farey tree, which is generated by the Farey map acting on the unit interval [0,1], or more 
precisely, on the point x = _l/2. It consists, at each level, of a subset of the Farey fractions. (For more details on 
these matters, see |3| and G3)- The Farey map is defined as 


/ fo{x) = x/(l - x) , if 0 < x < 1/2, 
\ fi( x ) = (1 — x)/x , if 1/2 < x < 1 . 


(23) 


We denote the inverses by F 0 (x) = / 0 *(x) = x/(l + x) and F\(x) = fi x (x) = 1/(1 + x). The associated Ruelle- 
Perron-Frobenius transfer operator is then formally given by (note the resemblance to (11(111 1 


Kp<Kx) = iFo'WlUiFoix^ + lFi'Wf^F^x)) 


1 + x 


1 + X 


(1 + x ) 2 / 3 l 

Therefore, the /c-fold iterated operator K, 1 ^ ip(x) consists of 2 k terms of the form 

I {F T1 oF T2 o...oF Tk )'(x)f<p{F Tl oF T2 o...oF Tk (x)) 


(24) 


(25) 


with Tj £ {0,1}. As we are dealing with iterations of Mobius transformations of the form with determinant ±1, 
we can alternatively consider multiplication of the associated matrices. We find for instance 


KyjiO) = ^(cx + d)^ = ^l(x)|Mi, 

{Tj} 


(26) 


2=1 

where c and d are just the bottom left and right entries, respectively, of the matrix 

k 


Mj = F T;i where F 0 = f | and F ± = f ^ j 


.7=1 


(27) 


Note that Aq = Fq and F\ = PA\. 
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When we apply K.p to the constant function l(ir) we obtain 2(1 + x)~ 2 ^. That is exactly twice the initial condition 
of the partition function Z k (x) (see @). In addition fCg increases the level k of the partition function Zk(x) by one 
as follows from ifTTill and d. Thus 

IC k 0 l(x) = 2Z k - 1 (x). (28) 

d is one of our main results. A connection of this type follows from and |3], but it is less direct, and not valid 
for all x values. Next, we consider some of the consequences of (|28|) . 

First, we note that for x = 0, d connects the Knauf model 0 and the transfer operator Kg: 

^l(x)U =0 = 2Zf_ 1 (2/3). (29) 

Now pj defines an operator C(2/l) whose non-degenerate leading eigenvalue A(/3) gives the free energy of the “grand 
canonical” partition function Zk(x = 1,(3) and the “canonical” case Zk(x = 0,(3) as in iBId) below. It also connects 
the largest eigenvalue of C(2(3) with the largest eigenvalue of the equation 

\{(3)f(x) = f{x + 1) + x~ 2f3 f{ 1 + 1/a;), (30) 

which is directly related to P . In fact the proof uses a Taylor series expansion of (fix) (in 11241) 1 at x = 1. However, 
the connection of the partition functions and d or d in 4| is not so direct. 

A connection between the spectrum of the operator C(2/3) and the spectrum of (1241) can be made, but the situation 
is complicated (in part because the spectrum determined in Q is on the space of functions of bounded variation), and 
not really germane to our purpose here, and so will be omitted. 

Next consider d for (3 > (3 C = 1. In that case, one has 

Km JC*l(aO| s=0 = Km 2 Z^(2(3) = 2 C( ^~ ) 1) , (31) 

where £ is the Riemann zeta-function (the second equality is shown in ]2|). This result has not appeared previously, 
to our knowledge. 

For (3 < (3 C = 1, the leading eigenvalue A(/3) > 1 of K,p [H is non-degenerate and belongs to the discrete spectrum. 
Since the corresponding eigenvector is of definite sign, it has a non-zero projection onto l(cc). Thus we can define 
aix , (3) as 


/ z kix,(3)^_ 

“ (x ' /3) = i” ~!?W 


(32) 


Note that since the spectrum of JCp is independent of x, the free energy 

-1 r ln Z ki x ,(3) -1 

f(P) : = J im - \ - = 

jj k—>oo K fj 

depends only on the inverse temperature (3 mu Thus (as we have already noted for x = 0 in [l5|l the phase 
transition is second-order for all x > 0. This follows from the result of Prellberg B©], 


(33) 


* [1 + 0 ( 1 )], ( 3 ^ r 


(34) 


ln(l - P) 

where c > 0 (for more discussion about the phase transition see 0 )- 

Now since all terms in Zkix,(3) are positive, and the matrix Aq + 1 is included (see 0 ), one has (recall that x > 0) 


(! + (*:+ l)x)~ 20 < Z k ix,(3) <Z k iO,(3)= Z? (2/3). 


(35) 


Thus, since the Knauf free energy vanishes for (3 > (3 C , so must the free energy obtained from Z k ix,(3). Furthermore, 
since the leading eigenvalue of ICp is A(/?) = 1 for all (3 > (3 C we can write for all temperatures 

fi x ,(3) =( 36 ) 

We have shown elsewhere [lfj that the free energy of the Knauf model, the Farey tree model and the Farey fraction 
spin chain of Kleban and Ozltik are the same for all temperatures and are also given by m- 
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Note that the leading eigenvalue changes its character at the critical point. Above the critical temperature it 
belongs to a discrete spectrum and below the critical temperature it is the upper limit of the continuous spectrum 
(for more details about the spectrum see |f|). 

The sub-leading eigenvalue in the spectrum is equal to one for P < p c . This is consistent with our previous results 
in based on scaling and renormalization group arguments. For a one-dimensional system the scaling arguments 
provide the relation between the singular part of free energy f s and correlation length £ 

fs oc i. (37) 

(The correlation length is essentially the distance over which the spin-spin correlation function varies, i. e. is not 
constant.) If we assume that our partition function goes as 

Z k (x, P) = A k a(x) + A^ai(a:) + ..., (38) 


we obtain, using i®, 


fs oc In A, 


and from the definition of the correlation length 


•e 


c 

ln(A/Ar) 


(39) 


(40) 


where C is a positive constant. This implies that the sub-leading eigenvalue Ai(/3) = 1 for P < /3 C , consistent with 
Prellberg’s results. 

In addition, note that from and the evenness of Z k , it follows that the eigenfunction a(x,/3) is even 

a(x, P) = x~ 2/3 a(l/x, P). (41) 


Using this fact and m we can write 


A (P)a(x, p) - a(x + 1, p) + (1 + x) 0 a ^ + ^ , PJ . 

(42) 

Note that, as remarked in 3 (and using the evenness of a(x,P)), 14211 is a generalization of the Lewis three-term 
equation, which has been extensively studied in number theory in the context of the Selberg £-function (@,0)' 
the Lewis case, solutions with A = 1 are of interest, and P may be complex. 

Applying 1411) and 1421) with x = 0 and x = 1 we obtain 

a (l, P) = (A(/?) - l)a(0, P), 

(43) 

and 


a(2,P) = ^a(l,P) = ^(A (p) - l)a(0,P), 

(44) 


respectively. We will make extensive use of (TTTtl and G3J below. 


V. EXPECTATION VALUES-PRELIMINARIES 


In this section we consider various spin expectation values for Knauf spin chains. (The remarks just below 121 1) 
define these quantities.) Making use of the spin flip behavior and recursion relations proved above, we obtain a few 
results, but our main purpose is to set the stage for the expectation value calculations of the next section. 

First, consider the expectation value for spin up 



l r 


Z l k U (x) 

Zk{x) 


(45) 
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(46) 


and similarly for spin down. By using (O, m and m we find 



Zi(x)\Aq(A$ + Ai) r 

Zi(x)\A 0 (A 0 + Ai) r + Zi(x)\Aq(Aq + Ai) r P 


We now relate the two terms in the denominator, at least for some values of x. We already know from lt22ll that for 
x = 1 these terms are equal. There is a simple explanation for this. Multiplying any matrix Mi by P on the right 
just exchanges its columns, and m is clearly invariant under exchange of columns for x = 1. Thus the probability 
to find spin up (or down) at any location on the spin chain with x = 1 is 

= (^> 1 ^Jx=i = ^ ' ( 47 ) 

l r l r 


Thus, in this case, due to the spin-flip symmetry, there are no finite size or edge effects at all (the result is valid for 
all l, r > 0). The situation is very different for x = 0, as we will see. 

Although our spin chains are defined in terms of matrices, one can also investigate their Hamiltonians. Generally, 
these are not very useful, since they include long-range many-body interactions between the spins (see (0) or l(nfl 
for definitions and explanations of these matters). However, it is known that for the “grand canonical” spin chain, all 
interactions are even and ferromagnetic (i. e. favoring aligned spins) Q. Therefore any expectation value involving 
an odd number of spins must vanish. Since x = 1 corresponds to this spin chain, GZJ is exactly what one expects. 

Now consider x = 0 (the Knauf model of (7J). The partition function at level k is 

Z k ( 0) = (Z,0r)|A o (A o + A{f + Z,(a;)|4)(4> + AtfP) | x=0 = 2Z‘ T ”(0) + Z t { 1) - Z t ( 0). (48) 

This result may be proven directly from the structure of the Farey fractions together with the action of the matrix 
P. However, we will show it by using 11 Mil and m . First, note that ZjJ° r ( 0) = Zi( 0) and Z l k ° r ( 0) = Zi( 1). Next, 
express m as 


Z‘ k ^(x) = (1 +x)-^Z^ (t^) +^ T _r (* + !)• (49) 

Now for x = 0 iRTill becomes 

Z^(0) = Z l (0)+ zf'-'-'Q.). (50) 

i=l -\-1 

Similarly we find 

k -1 

Z^(0) = Z l (l)+ (51) 

i=l -\-1 


Adding the above expressions (seeJ2J) and using the fact that Z\} r { 1) = Z‘ f } r ( 1) (see iF^ll l leads to (HH1) . 
By making use of iTTfTl) and the expectation value for spin up at x = 0 can be written as 


where 



1 

2 - I< ’ 


Zi(0)-Zt(T) 

Zi(x)\Aq(Aq + Ai) r \ x=0 


(52) 


(53) 


Now Zi(x)\A 0 (A 0 + Ai) r |a, =0 > Zi( 0) > Zi( 1) > 0 for all l,r > 0 and f3 > 0. The first inequality follows immediately 
from the fact that the sum Zi(x)\A 0 (A 0 + Ai) r \ x= o of positive terms includes the term Zi(x)\A t 0 +1 \ x= q = Zi( 0). The 
second inequality follows directly from the monotonicity (in x) of Zi(x). Therefore 0 < K < 1, where I\ = 0 can 
occur if the denominator of 15311 diverges. (This happens when r —> oo, see 1661 below.) Thus the spin at any position 
for temperature T < oo has, in general, greater probability to be up than down 





0j 


l 


l 


(54) 






where it should be realized that equality only holds in the special case K = 0. 

In the “normal” situation, i. e. when equality does not hold, im may be regarded as an effect of the “hidden” spin 
up on the left, i.e. the initial condition (1 + x)~ 2 ^ = l(ir)|Ao, which breaks spin-flip symmetry. A slightly different 
point of view involves the spin interactions. For x = 0, i. e. the (canonical) Knauf model, these are all ferromagnetic 
and include terms with an odd number of spins 0 ii. isj shows that the odd interactons can be sufficient to favor 
an up spin. The interactions also give rise to some rather subtle effects in certain other expectation values, as we will 
see below. 

Now we consider the two-spin correlation function. Let 



Z l J n (x) | Ao(Ao + A\) r 
Zk+r+l {x) 


(55) 


where as before k = l + n + 1. 

The partition function Zk+ r +i(x) for a spin chain of length l + n+r+2 can be divided into four terms (corresponding 
to the four possible configurations of two spins) 

Z‘ k u (x)\A(A 0 + Atf 

and 

Z{x)\Ai(AQ + Ai) r 

where i £ {0,1}. Using the matrix P (see (J22J) gives 


Z l t } n {x)\Ai[AQ + Ai) r — Z (a:)|^4i+i(mod2)(^4o + Ai) r P. 


(56) 


Now m shows that for x = 1 each spin has equal probability to be up or down without any edge or finite size effects 
(i.e. for any l, r £ Zj). Thus we can expect that e.g. the expectation value for two spins up is the same as for two 
spins down. In fact jnni and (I5H give immediately 


and 




X = 1 


l 


n r 



(57) 


(58) 


where l, n, r £ 

In the case of one spin m shows that the expectation value does not change under translation of the spin. The 
two spin expectation value is not translationally invariant but it does have the following symmetry 

( 59 ) 

l n r r n l 

This follows on rewriting the 1. h. s. of as 

ZjJ n (x)\Ai(Ao + Ai) r = (1 + x) ~ 3 \(Aq + Ai) 1 Aq(Aq + Ai) n A\(A q + Ai) r 

<2^ l + n, + r 

= ^ [(a + c)x + b + d]~ 213 , (60) 

i=i 

where a, b , c, d are entries of the ith matrix from the set {Aq + Ai) 1 Aq(Aq + Ai) n Ai(Ao + Ai) r . Thus for x = 1 the 
sum does not change under matrix transposition and we get This result also follows from the proof that the 

interactions in the “grand canonical” spin chain also have the symmetry (1591) (Lemma 4.8 in [§]). 
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VI. EXPECTATION VALUES-RESULTS 


In this section we calculate several spin expectation values for Knauf spin chains. These results are all new, and 
are the first calculations of such quantities, to our knowledge. We find that they are expressed as simple functions of 
the free energy / (or correlation length £). 

The methods used in the previous section (and this one as well!) are only of use when one can come up with a finite, 
closed set of equations. For the expectation value of a general set of spins, this is not the case. Since the matrices 
representing the spins operate “from the right” most of our results are expectation values involving a finite number 
of spins fixed at or at a finite distance from the right hand side of the spin chain. 

First note that by El, at x = 1, the spin expectation value has no edge or finite size effects. Thus, allowing l —> oo 
and r —> oo, a spin up (down) still has probability one half. By contrast, 15211 shows that there may be such effects 
for x = 0. This indeed occurs, as we now proceed to demonstrate. 

In order to see the edge effect at the right side of an infinitely long chain we go back to H51) and let / —*• oo . Using 
(|32|) and the properties of the eigenfunction a(x) we get 

A r+1 a(0) = 2 (i(x)\Aq(Aq + A{) r \ x= o + a(l) — a(0) (61) 

for all r > 0. Thus, we can write the expectation value for a spin r + 1 from the right of the infinitely long chain using 

61 




(62) 


where the eigenvalue A(/3) € (1,2] for P G [0, /5 C ). Note that a similar expression for the spin down expectation value 
follows since their sum must be one. 

Now recall (see m) that A is given directly in terms of the free energy via A = e The free energy is a 
non-increasing function of the temperature for (3 < (3 C . Hence, for any fixed r, the expectation value m decreases 
monotonically to 1/2 as T —> oo. In fact, all our results are consistent with a product distribution in this limit, 
i. e. the probability of a given spin being up or down is 1/2. Note also that A may also be expressed in terms of the 
correlation length £ (see HOI and recall that Ai = 1 for (3 < f3 c ). 

From a physical point of view, it is also interesting to compare m and El, which, as mentioned, may be attributed 
to long-range interactions between an odd number of spins. |62|) shows that their effects are felt even infinitely far 
from the initial (“hidden”) up spin. This is particularly interesting, since 0 proves that even though there are 
odd (ferromagnetic) interactions at x = 0, any individual interaction term vanishes in the limit of an infinitely long 
chain. Thus El shows that certain cumulative effects of the odd interactions remain in this limit, even though each 
individual interaction goes to zero. 

It is also of note that El, as well as various expressions that we will derive shortly, give expectation values as 
simple polynomials in A, which is itself exponential in the free energy / or correlation length £, as mentioned. 

We can use El at the critical temperature (where A ((3 C ) = 1) by taking the limit (3 —> f3 c . Then the probability of 
a spin up is 1 for any finite distance r from the right (this can also be shown directly from 152H since K —> 1 when 
(3 —*• p c and then l — > oo). On the other hand for any (3 < p c the spin up (or down!) probability goes to one half as 
r —> oo. 

Note that El also gives the right edge correlation length as 


£r = 


1 

In A 


1 

Ts' 


This equation directly relates edge and bulk behavior. Since the bulk correlation length £ oc -j- (see (T6])), 


(63) 


oc £ oc 


lnt 


(64) 


as [3 —> Pc, where e = % — 1. 

Now consider the limit r —> oo, keeping l finite. Using El we can write 


lim 

r—► oo 


Zi(x)\Aq(Aq + Ai) r \ x =o 

W 



(65) 
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for any A > 1 (i.e. [3 < j3 c ). From ( 17911 below we see that 0 < a(0) < oo for A G (1, 2]. Using 15211 and © we then 
obtain 

T v^)*=o = \ (66) 

l OO 


for all Z > 0 and (3 < (3 C . Thus the left edge effects on one spin vanish. From a physical point of view, this is quite 
interesting. The “hidden” spin up on the left, or equivalently the long-range odd ferromagnetic interactions @ have 
no effect in an infinite chain when the spin in question is only a finite distance from the “hidden” spin. By contrast, 
when it is infinitely far away but at a finite distance from the right edge, there is an effect (see (1621) 1. However, we will 
see that this effect is removed if one fixes spins on the right hand end of the chain in certain specific configurations. 

Next we consider the two spin correlation function, in the limit where the left part of the spin chain goes to infinity. 
Using (O, <G2> and (1531) we find 



a(x)\A 0 (A 0 + Hi) "Ho (Ho + Hi) r 
A n+r+2 a(x) 


(67) 


It is convenient to define two functions of x and f3 (where, as for a(x), we do not explicitly indicate the (3 dependence), 



U n {x) = a(a:)|H 0 (Ho + Hi) n 

(68) 

for spin up and similarly for spin down 

D n {x) = aOr)|Hi(H 0 + Hi)”. 

(69) 

Clearly for all n > 0 and 0 < (3 < (3 C = 1 

U n {x) + Dnix) = A” +1 a(x). 

(70) 

Using {HI and (^J 

Unix) = x~ 2l3 D„il/x). 

(71) 

Now return for a moment to the one-spin 

expectation value. We can write 


< s 

1 \ U n ix) 

^ ^ Unix) + Dnix)' 

(72) 


For x = 1 it immediately follows from ra that 

= (73) 

oo n 


as already shown by llT7ll . Note also that (pll. for one spin at x 


0, follows from equation m which we rewrite as 


Un(0)= Q(A" +1 -A) + l)a(0), 

and similarly 

D n (0)= Q(A” +1 +A)-l)a(0). 
Now return to equation (1571) for r —> oo. First we rewrite it as 



oo nr 


Un(x)\A 0 (A 0 + Ai) r 
A n+r+2 a{x) 

EUn&jcx + d)-^ 
A n+r+2 a(x) 


(74) 


(75) 


(76) 
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where the sum has 2 r terms, and a, 6, c and d are from AqM v = |“Jj. Note that we start with the matrix Aq and 
thus < 1 for all x £ R+. Now 


U n (x) 



(ex + d)j 2/j 
\{(3) k 


(77) 


where j £ {A 0 (A 0 + A\) k A 0 (A 0 + Ai) n }. Since Cj, di > 0 for all i, it follows that U n (x) is non-increasing with 
x. Thus we can write 


U n (l)Z r (x) 

A n+r+2 a(x) 



U n (0)Z r (x) 

A n+r + 2 a(x) 


(78) 


for all r > 0 and all x £ Rj}. In the limit r —> oo we get, using m and U n (l) 


A n+1 a(l)/2 (see d and CD) 


(A-l) 


a(0, A) 
2A 



oo n oo 



2-A 
A" +1 


a(0,A) 

2A 


(79) 


Physically, since the correlation length £ = the n-dependence of the upper bound in d is what one expects for 
the correlation function itself. We have not been able to prove this, however. 

Now we calculate some results at the right hand edge, i.e. for finite r, with x = 0 (the “canonical” case). When 
r = 0 we have from CD and m 


Similarly 


as well as (see lpl>) 


and 



oo n 


Un{ 0) _/ 2-A\ 1 

A ra+2 a(0) ^ A ” +1 ) 2A 



oo n 


D n ( 0) / 2 — A\ 1 

A ra+2 a(0) ^ A n+1 ) 2A’ 



oo n 


Un{ 1) = A-l 
A” +2 a(0) 2A 



oo n 


Dn{ 1) = A-l 

A” +2 a(0) 2A 


(80) 


(81) 


(82) 


(83) 


It is easy to see, for instance, that the sum of and (El is the same as El, and that El and El sum to El 
with r = 0. Similar checks verify other sums of the four equations just above. Further, since the interactions are 
ferromagnetic, El must be the largest of the three, and this is easily verified as well (recall that 1 < A < 2). 

It is interesting that both El and El are completely independent of the spin separation n, and equal to each 
other for (3 < (3 C . Thus a down-spin at the right hand edge completely cancels the lack of spin symmetry seen in El , 
and does so for all (3 < (3 C . (In fact this holds for f3 > (3 C as well, since all spins are up in an infinite chain, so both 
El and El vanish.) We comment further on this after deriving some more general results. 

The results in the paragraph above are based on our knowledge of U n (x) and D n (x) at the two values x = 0 and 
x = 1. It is easy to find generalizations. We need combinations of spins for which the corresponding product of 
matrices A 0 and A\ has b = 0, so that x = 0 is preserved, or b = 1 and d = 1, so that x = 0 maps to x = 1. This 
is true for chain of Aq matrices of any length and chains starting with A\ following by a chain of Aq matrices of any 
length. These two cases give us certain expectation values with r spins on the rhs fixed: 

'O T •••]••■ t >*=o = { 1+ W+r) 2W (84) 


OO 
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and 


(85) 


1T... T... T)*=o = 


A- 1 
2A r 


Similarly we get 


... I ... t...T-.-T>*=0 = 1- 


oo n 


2 - A\ 1 


A n+1 J 2A r 


( 86 ) 


and 




oo n 


A- 1 
2A r ' 


(87) 


Note that E3 and IH7I generalize eh and (El in that they are both independent of n and equal to each other. 
Thus the restoration of spin symmetry already seen, which holds at any temperature and for any separation n, is also 
valid for any r. This seems very curious and nontrivial. Some understanding can be gained by considering recent 
ideas about the mechanism underlying the phase transition E3 According to this work, the transition is due to the 
condensation of clusters of spins of exactly the type on the right hand edge in El and El- This is consistent with 
our results here, since it shows that such clusters restore the spin symmetry which is broken by the “hidden” up spin 
on the left hand edge of the chain, at least for the particular expectation values investigated. 

Note that to calculate any of El - El for x = 1 would require knowing the four values J7 n (l/2), U n ( 2), D n { 1/2) 
and D n ( 2). Using (T7T1) accounts for two of these, in addition d removes one more, but one is left with one unknown 
value. For general x 7 one has four unknown quantities. 

Finally, summing El and El or (1531) and El gives rise to, respectively, 


• • T •••?•■• T)x=o = 


( 88 ) 


and 


I T ■ ■ • T ■ • • t)x=o = 


A- 1 


(89) 


As one approaches the phase transition, A —> 1, so that El goes to 1 while El approaches 0. This suggests that in 
an infinite chain exactly at the transition, the only state with non-zero probability has all spins up. This would not 
be surprising, since the same property holds below the transition (for /3 > /3 C ). 


VII. CONCLUSIONS 


In this paper we have extended our understanding of the statistical mechanical behavior of the Farev spin chains. 
Our main tool is a generalization of the “number theoretic” partition function studied by Knauf { 2 , 13, 4 L| • 

By introducing a new parameter, we are able to derive recursion relations on the length of the spin chain (or 
equivalently, the level of the Farey fractions). These relations are generalizations of the Lewis three-term equation of 
number theory S0 Using them and the behavior of the system under spin-flip transformations, we find new results. 
In particular, we prove a new and simple connection between the Ruelle-Perron-Frobenius transfer operator studied by 
Prellberg [§, !J, ID] in a model of dynamical systems with intermittency and our generalized partition function. This 
connection implies that all our models have the same free energy (and hence the same thermodynamics), independent 
of the value of x. It also implies that the correlation length satisfies the prediction of scaling theory. 

In addition, we are able to calculate certain spin expectation values for both finite and infinite spin chains. These 
results are the first such calculations, to our knowledge. In particular, we find that the expectation value of certain spin 
clusters and an independent spin at arbitrary distance is independent of the direction of the spin, at all temperatures 
above the transition. This holds even though the expectation value of the the independent spin by itself does depend 
on direction. Thus, in this sense, the spin cluster removes the spin asymmetry of the system. This behavior appears 
to be related to the mechanism of the phase transition. 
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